Spiral Exchange Coupling in Trilayer Magnetic Junction Mediated by 
Diluted-Magnetic-Semiconductor Thin Film 
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We revisit the non-collinear exchange coupling across the trilayer magnetic junction mediated by 
the diluted-magnetic-semiconductor thin film. By numerical approaches, we investigate the spiral 
angle between the ferromagnetic layers extensively in the parameter space. In contrast to previous 
study, we discovered the important role of spin relaxation, which tends to favor spiral exchange over 
the oscillatory Ruderman-Kittel-Kasuya-Yosida interaction. Finally, we discuss the physics origins 
of these two types of magnetic interactions. 

PACS numbers: 75.50.Pp, 75.30.Et, 72.25.Rb, 75.70.Cn 
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The goal to merge the functionalities of informa- 
tion processing and data storage in one single mate- 
rial has charmed researchers in the field of spintronics 
for many years. Q One of the promising candidates 
is the diluted magnetic semiconductor (DMS), made of 
the III-V host semiconductor doped with transition met- 
als, such as (Ga,Mn)As. The ferromagnetic order in 
DMS0, 0, 0, |a, arises from the aligned magnetic mo- 
ments of the transition metals, mediated by itinerant car- 
riers in the semiconducting bands and thus can be easily 
manipulated by electrical means. This unique feature has 
attracted enormous interests both in academic research 
and potential industrial applications. 

In addition to its potential applications for making 
the next-generation transistors in spintronics, DMS also 
brings up surprises in the more conventional magnetic 
multilayers as shown in Fig. ^ One of the authors stud- 
ied the F/DMS/F trilayer magnetic junction, within the 
linear response theor y [al and the self-consistent Green's 
function approach |jl lift ITlj. and found interesting spi- 
ral exchange coupling between the ferromagnets. Their 
numerical results show that the spiral exchange always 
beats the Ruderman-Kittel-Kasuya-Yosida (RKKY) in- 
teraction, rendering it into decorative ripples on the spi- 
ral backbones. Being puzzled by the dominance of the 
spiral exchange, we revisited the problem and found the 
important ingredient overlooked in the previous study - 
the proper inclusion of the spin relaxation rate. 

In this Letter, employing the same model described in 
Ref.Q, we studied the trilayer magnetic junction within 
linear response theory. We carried out extensive numeri- 
cal computations in different parameter regimes, in par- 
ticular, with different spin relaxation rates r) (which was 
fixed in the previous study). We found that the compe- 
tition between the spiral exchange and the RKKY inter- 
action sensitively depends on the strength of spin relax- 
ation. In the extreme ballistic limit, the RKKY interac- 
tion prevails, while the spiral exchange starts to take over 
the leading role when approaching the diffusive regime. 
The value of 77 used in Ref. 8| happened to set the sys- 
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FIG. 1: Schematic figure for the F/DMS/F trilayer magnetic 
junction, where the spiral angle 3 (r) denotes the relative an- 
gle between the two ferromagnets. 



tern in the diffusive regime (see below) where the spiral 
exchange dominates. In general, depending on the spin 
relaxation rate, it is possible to observe both types of 
magnetic behaviors! 

Before diving into numerical details, a simple argument 
from the single-particle picture would help readers to un- 
derstand the origin of the spiral exchange. In fact, the ar- 
gument closely parallels Datta-Das' original proposal [l^ 
for a spin field-effect transistor. Imagine an electron with 
+x spin orientation is injected into the DMS thin layer 
from the ferromagnet on the left-hand side. The two com- 
ponents of the wave function evolve differently because 
of the Zeeman splitting in DMS, 



*(0) 



V2 V 1 



*(r) 



x/2 (c 



(1) 



here r is the thickness of the DMS thin layer. The phase 
difference between the spinor components indicates that 
the carrier-mediated exchange coupling is non-collinear 
and tends to align the other ferromagnet on the right 
at a different angle, 9 s (r) = —2irr/A s = — (kp^ — kpi)r. 
However, this simple single-particle picture may not be 
the whole story because it does not capture the quantum 
interferences from different patches of the whole Fermi 
surface, which give rise to the oscillatory RKKY interac- 
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tion. Therefore, to pin down the dominant magnetic in- 
teraction, one needs to resort to the more formal Green's 
function approach. 

After integrating out the itinerant carriers, it can be 
shown that the effective exchange coupling between the 
ferromagnets is proportional to the static spin suscepti- 
bility, 



X ij (r) 



^{{i[a\r, t),o*(p,0)]))dt, (2) 



where a 1 is the spin density in DMS layer and the double 
bracket denotes thermal and quantum mechanical aver- 
ages. Note that the vector dependence can be dropped 
because of the external SO (3) symmetry for the spatial 
coordinates (not to be confused with the internal SO (2) 
symmetry for the spinor). Furthermore, a phenomenol- 
ogy parameter r\ is introduced to describe the spin relax- 
ation in DMS layer. 

In the presence of the finite Zeeman gap, the spinor has 
an internal SO (2) symmetry. Consider a 7r/2-rotation 
along the z-axis (always chosen to be the quantization 
axis for the Zeeman splitting). It changes (S x ,S y ) — > 
(S y ,—S x ) and leaves S z intact, implying \ xy = — X yx 
and x xx = X vy ■ Similar argument leads to x xz = = 
X yz ■ For general geometries, the orientation of the pinned 
ferromagnet can be — (sin 6* cos (/>, sin sin 0, cos 9). 
Mediated by the itinerant carriers, it will lock the free 
magnetic moment on the right-hand side in the direction 
of n R , 
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Since our focus here is to determine the dominance of 
different magnetic interactions, we choose the simplest 
geometry in Fig. ^ by setting = (1,0,0). Thus, it 
is clear that tir can be described by the spiral angle 
6> s (r) = -tan _1 [x xy (r)/x xx (r)]. To find out the spiral 
angle, we use the well-known trick to compute the com- 
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FIG. 2: Spiral angle s (r) for (a) k A /k F = 0.5 (b) k A /k F = 
1.0 (c) k A /k F — 1.3 in the ballistic regime with k v /k F = 0.23. 



obvious time-reversal symmetry within each spin flavor. 



plex susceptibility x = X xx + i-X 1 
desired real and imaginary parts, 



then extract the ^ gives rise to the modified Onsager relation \ 
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where = k 2 /2m* is the dispersion for itinerant carriers 
and A is the Zeeman gap. The Fermi-Dirac functions for 
the itinerant carriers are f a {tk) — [exp(e^ — a A — /i) + 
1] j with a = ± corresponding to up/down spin flavors. 

When carrying out numerical calculations for the spin 
susceptibility, we found that Xxy ~ V 2 vanishes as r\ — > 
+ . It implies that RKKY dominates in the ideal ballistic 
limit. Actually, in this limit, the spin-spin commutator 
in Eq. Q can be computed analytically and is indeed 
zero. We later realized that this result is due to a less 



ideal ballistic limit. In Ref.|8|], the authors only chose 
one particular rj for their numerical computations (which 
happens to be in the diffusive regime). Therefore, their 
conclusions are incorrect in the extreme ballistic limit. 

In realistic materials, the spin relaxation rate rj is fi- 
nite and the symmetry constraint no longer applies. By 
gradually increasing rj from the ballistic to the diffusive 
regimes, we show how the trend of the spiral angle 9 s (r) 
changes with it, as summarized in Figs. [21 and 03 For 
convenience, we introduce &a = V2m*A to denote the 
inverse-length scale for the Zeeman gap and k n = y/2~m*rj 
for spin relaxation. The total density of the itinerant car- 
riers is also converted to hp. In DMS, &a has a sensitive 
dependence on the temperature through the Zeeman gap 
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FIG. 3: 6s(r) for k v /kF — 1.0 in the diffusion regime: (a) 
(WM = 0.5, (b) (k A /k F ) = 1.0. 

and can be as large as (or larger than) kp (roughly 1/nm 
in typical DMS materials) at low temperatures. Further- 
more, the spin lifetime is about 10 -2 picosecond, giving 
the ratio k v /kp w 0.23. 

Making use of this rough estimate, we compute the 
spiral angle with different Zeeman splitting )ca , shown in 
Fig. When the polarization is small {k^/kp = 0.5), 
the signature of RKKY oscillations is rather obvious. 
On the other hand, when the carriers are fully polarized 
{kA/kF = 1.3), the spiral rotation becomes transparent 
and the RKKY oscillations, although visible, are sup- 
pressed into minor ripples. Careful analysis shows that 
these complicated patterns for the spiral angle can be 
characterized by two length scales, 

ARKKY — 1 —j , *s - 1 ; ■ (O) 

Kfi + kp[ Kfi — Kfi 

These length scales originate from the low-energy spin ex- 
citations. Consider the particle- hole excitation by kicking 
a spin-down electron to the spin-up band. It will carry 
momentum p = kp-\ — kp^. The length scale A s arises 
from excitations with kp^ an d fcpj. parallel, while Xrkky 
comes from those with antiparallel momenta. The com- 
petition between these two type of spin excitations leads 



to the complicated patterns of the spiral angle s {r). 

To explore the role of spin relaxation, we also extend 
our calculations to the diffusive regime k v /kp = 1 as 
shown in Fig. The spiral exchange becomes signifi- 
cantly enhanced, rendering the oscillatory parts into dec- 
orative ripples. This can be understood as the decease of 
the time-reversal symmetry which is broken by the large 
spin relaxation rate. Note that the RKKY oscillations 
come from quantum interferences between the patches 
of the Fermi surfaces related by time-reversal symmetry. 
Thus, by breaking the symmetry, it is very efficient to 
weaken the RKKY interaction, as demonstrated by our 
numerical results. We would like to emphasize that the 
previous study fails to recognize the important role of the 
spin relaxation and thus misses out the subtle competi- 
tion between the spiral exchange and the RKKY inter- 
action completely. 

In conclusion, we demonstrate the important role of 
spin relaxation in the trilayer magnetic junction and com- 
pute the spiral angle between the ferromagnetic layers 
with different carrier concentrations, (temperature de- 
pendent) magnetizations and spin relaxation rates. Our 
numerical studies show the non-collinear coupling across 
a DMS thin film is important and will play a crucial 
role for magnetic junctions at nanoscale. Finally, grant 
supports from National Science Council in Taiwan are 
greatly appreciated. 
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